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Abstract. We study the phase ordering colloids suspended in a thermotropic nematic liquid crystal below 
the clearing point T n i and the resulting aggregated structure. Small (150 nm) PMMA particles are dispersed 
in a classical liquid crystal matrix, 5CB or MBBA. With the help of confocal microscopy we show that small 
colloid particles densely aggregate on thin interfaces surrounding large volumes of clean nematic liquid, 
thus forming an open cellular structure, with the characteristic size of 10 — 100 /im inversely proportional 
to the colloid concentration. A simple theoretical model, based on the Landau mean-field treatment, is 
developed to describe the continuous phase separation and the mechanism of cellular structure formation. 

PACS. 61.30.-v Liquid crystals. - 82.70.-y Disperse systems. - 64.74.+g Solubility, segregation and mixing; 
phase separation. 

1 Introduction 

The phase separation and ordering of multi-phase systems, 
and their resulting physical properties, have been the sub- 
ject of much active research over the past twenty years. 
Such a high interest in this field arises from a combina- 
tion of challenging, fundamental physics (of phase equi- 
libria and dynamics) with a large number of viable ap- 
plications. Research in liquid crystals has also remained 
active for decades, for similar fundamental and practical 
reasons. However, the applications of thermotropic liquid 

. crystals have always focussed on display technology, some- 
what overlooking their unique mechanical properties. In 
the early 1990's a new field, now called "liquid crystal 
colloids" , was opened to research by an experimental and 
theoretical study of suspensions of colloidal particles in a 

. lyotropic liquid crystal by Poulin et al. Q. The key idea 

1 of introducing a large amount of mobile interface into the 
liquid crystal, which generates nontrivial topological con- 
straints and singularities, has been fruitfully explored by 
further experimental and theoretical work 

The problem studied is that of a spherical particle with 
nematic director anchored on its surface, suspended in a 
nematic liquid crystal matrix. When the director is an- 
chored rigidly, a closed inner surface creates a topological 
mismatch between the director field n(r) on the parti- 
cle surface and the uniform director at large distances 
from it. This mismatch leads to topological defects, i.e. 
regions where the liquid crystal order and the continuity 
of n break down. A connected closed surface represents 
a point topological charge N = 1. Since the overall sam- 
ple has to be topologically neutral, the charge produced 



by the inner surface must be balanced by an associated 
opposite charge. There are two basic possibilities. The as- 
sumption, that a spherical particle in a quadrupolar ne- 
matic medium creates deformations that are quadrupolar 
too, results in the mismatch balanced by a circular loop 
of disclination with linear strength N = —1/2 and overall 
point charge N = — 1 , in the equatorial plane of the parti- 
cle: a "Saturn ring" M, see fig]l|(a). The other possibility 
is a dipolar configuration with a satellite point defect, a 
monopole with N = — 1 near one of the particle poles ||, 
see figj^(b). It appears that, although both situations arc 
possible, the dipolar structure with satellite defect occurs 
more readily, e.g. in experiments aM for water droplets 
dispersed in a thermotropic nematic matrix. When sev- 
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Fig. 1. The calculated director field around a single spheri- 
cal particle in a uniform nematic matrix. For rigid anchoring 
there are two possibilities: (a) the quadrupolar disclination ring 
structure, or (b) the dipolar structure with a satellite monopole 
defect. In the case of weak anchoring (c) the director field has 
the quadrupolar symmetry but does not have topological sin- 
gularities. 
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eral such droplets come to a near-contact, the effect of 
chaining occurs, with satellite defects sandwiched between 
inclusion particles and shared between them, providing a 
strong and anisotropic interaction force. 

On the other hand, when the strength of director an- 
choring is rather weak, the particle introduces only a small 
perturbation into the nematic matrix, see fig.[l](c). Albeit 
small, the non-uniform distortion of n is spreading over 
long distances from the particle - the exact solution for 
the angle 9 of director deviation from the uniform axis n Q 
can be obtained in terms of the multipole expansion and 
reads, e.g. H (in spherical polar coordinates), 
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In this expression, R is the particle radius and the nematic 
constants W and K are the anchoring energy on its sur- 
face and the Frank elastic constant respectively ||. Their 
typical values in a thermotropic nematic liquid crystal are 
K ~ 10 _11 J/m and W ~ 10~ 6 J/m 2 , for homeotropic di- 
rector anchoring as in fig.[j]. 

The parameters in eq. (|lj) identify the important di- 
mensionless ratio which controls the outcome of the sin- 
gle particle behaviour. (WR/K) measures the relative ef- 
fect of the particle surface (characteristic energy ~ WR 2 ) 
and the bulk director deformations (Frank energy scale 
~ KR), and can be large or small. When WR/K <C 1 
the anchoring should be considered weak and not capa- 
ble of producing large deformations in the surrounding 
nematic matrix. In contrast, when WR/K 3> 1 rigid an- 
choring can be assumed and topological singularities re- 
sult. For typical parameters, particles of R ~ 5 — 10 fim 
represent the border between the two topologically dis- 
tinct regimes in fig.[l]: submicron colloid particles create 
small long range director distortions of quadrupolar sym- 
metry, eq. ([!]), while larger objects most frequently have a 
dipolar configuration with a satellite defect. 

The second physical idea relevant to our work is using 
the process of phase separation in a liquid crystal phase as 
a means of creating internal interfaces. The orientational 
symmetry breaking and the additional curvature elastic 
energy in the liquid crystal microconfined by such inter- 
faces should have a profound effect on the whole process 
of phase ordering. Both phase equilibrium and kinetics 
are altered in non-trivial ways due to the underlying frus- 
trated liquid crystalline order. 

Although the progress in the single-particle descrip- 
tion of liquid crystal colloids has been noticeable, little 
experimental work exists in this area of colloid collective 
behaviour. One of the main questions here is the macro- 
scopic morphology of a colloid and the change in rheo- 
logical behaviour on increasing concentration of particles, 
when the suspending matrix would possess a dense net- 
work of interaction forces creating high barriers to defor- 
mation. One may expect a glass-like freezing of motion at 
sufficiently high concentrations - an interesting possibility 
that needs to be thoroughly explored. 

The work of Poulin et al. on water droplets in a ne- 
matic matrix confirmed the ideas about topological defects 



and their long-range interactions - but they only studied 
the structure around a single or few very large droplets. In 
contrast, Tanaka et al. [^) studied phase ordering in mix- 
ture of a nematic 5CB and a surfactant, which forms very 
small nano-micelles. They reported some provocative re- 
sults on a phase they named "transparent nematic" . One 
could speculate that in this case the colloid parameter 
WR/K is so small that the mean-field description of ne- 
matic order would fail altogether. Recently, Zapotocky et 
al. [jl0| examined the colloid rheology of a cholesteric liq- 
uid crystal filled with silica particles of size ~ 1 /im. They 
have found a variety of rich behaviour due to topological 
defect-driven enhancement of dynamic elastic properties 
and the formation of a weak gel with the estimated static 
modulus of G' ~ 0.01 Pa. However, in their work small 
silica particles have aggregated into much bigger objects 
(floes) that produced strong topological forces and ended 
up in the nodes of a network of disclination lines, which 
has been carrying the clastic function of the resulting ma- 
terial. 

A simpler system has been discovered recently by Meeker 
et al. jn]. A "classical" liquid crystal colloid has been 
prepared by mixing a well-characterised thermotropic ne- 
matic liquid crystal, 5CB, with sterically-stabilised PMMA 
particles, imposing radial boundary conditions in the ne- 
matic matrix. Mixtures of 5—10 % particle volume fraction 
resulted in a soft solid with significant storage modulus, 
G ~ 10 3 — 10 5 Pa. Such a remarkable mechanical transfor- 
mation requires further investigation, which is reported in 
the companion paper ]T2| . Bright field microscope observa- 
tions by Meeker et al. suggest that the soft solid comprises 
a network of particle aggregates, formed by the exclusion 
of particles from emergent nematic domains as the mix- 
ture is cooled below the isotropic- nematic transition T n ;. 
However, their optical observations were limited to low 
particle volume fractions (< 5%), and depended on re- 
heating the samples back into the isotropic phase to view 
the particle aggregate structure, since direct optical study 
in the bulk of a highly non-uniform birefringent system is 
difficult due to the strong scattering of light: the nematic 
colloids are opaque below T nl . In this paper we use a confo- 
cal microscopy technique to study directly the structures 
in the bulk of the nematic liquid crystal colloid samples, 
for a wide range of particle concentrations (3 — 15 %). 

In addition to the remarkable mechanical response, 
the second aspect of interest in this system is the time 
required for the formation of an aggregated state below 
T n j. The characteristic time required for particle move- 
ment under the influence of long range attraction forces 
has been estimated in || and gives, for our example of 
particle size and concentration in the limit of weak an- 
choring, t - 10 3 ri(K/W 2 ) ~ 10 3 s (here rj ~ 0.1 Pa.s 
is the relevant viscosity coefficient of the nematic liquid 
crystal) . We shall see that a rigid gel is formed in a much 
shorter time after the colloid mixture is brought below the 
nematic transition point T n ;. It appears that the phenom- 
ena observed in our liquid crystal colloid system cannot be 
purely due to the topological defects and their networks: 
their elastic energy is not enough to explain the modu- 
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lus and, more importantly, there should be no topological 
defects at small WR/K. 

Here we study the formation process of the resulting 
near-solid liquid crystal colloid aggregates in more detail. 
We show that a process of continuous phase separation of 
small particles takes place in our system. This separation 
occurs immediately after the nematic transition and re- 
sults in the formation of a metastable but very long-lived 
rigid cellular structure of very thin interfaces, where the 
particles are densely packed, encapsulating large volumes 
of nearly pure nematic liquid. We develop a mean-field 
theoretical model describing the phase stability and the 
mechanism of cellular solid formation. 



2 Sample preparation and experimental 
methods 

In order to investigate the robustness of the process dis- 
covered by Meeker et al. jy], we study two materials: 
5CB and MBBA, both archetypical thermotropic nemat- 
ics |l3). The case of 4'-pentyl-4-cyanobiphenyl, abbrevi- 
ated as 5CB, obtained from Aldrich Chemicals Co., has 
already been examined in jyj. This pure material has a 
stable nematic phase below T n i = 35. 8C, followed by a 
crystal phase at T x = 22C. No crystallisation was observed 
down to 5C when the colloid particles were added. 

In contrast, the MBBA (Af-4-methoxybenzylidene-4'- 
butylaniline, from Aldrich), although one of the most com- 
mon nematic materials in the literature, has a practi- 
cal complication of being susceptible to hydrolysis. The 
C = N bond in the rigid molecular section dissociates in 
the presence of water. The reaction is reversible, mean- 
ing that at a given temperature and humidity there is an 
equilibrium balance of proper MBBA compound and its 
hydrolysis products, which act as impurities and weaken 
the nematic order. The effect is most noticeable in the 
depression of the phase transition points: in the "dry" 
MBBA the nematic transition point is at T$ = 47C fd| , 
followed by a crystal phase at T x — 22C, while we have 
used MBBA in hydrolysis equilibrium at ambient condi- 
tions, with T n i — 37. 3C and no crystallisation down to 
5C. We employed this as a model for a "dirty", weakly 
nematic liquid. 

Theoretical arguments suggest that the best phase sep- 
aration conditions are achieved for small colloid particles 
(the ones with WR/K <C 1, but not too small so that the 
reversal of this relation can occur near T n i), see section [| 
below. We use monodisperse PMMA spheres (polymethyl- 
metacrylate, bulk glass transition at T g = 105C) of radius 
R = 150 nm, with a polydispersity of about 0.04 (as de- 
termined by dynamic light scattering). Particles were cov- 
ered with chemically grafted poly-12-hydroxystearic acid 
(PHSA) chains (prepared by Dr. A. Schofield in Edin- 
burgh). On grafting, the short stabilising chains adopt a 
conformation radially extending from the grafting surface, 
making a "hairy" particle. In isotropic suspending liquids, 
these particles behave like almost ideal hard spheres. In 
addition, in the nematic suspending matrix, the grafted 



chains provide a homeotropic (radial) director anchoring 
with a typical strength W ~ 10~ 6 J/m . With the typical 
Frank constant K ~ 10 -11 J/m this makes the dimension- 
less colloid parameter to take the value WR/K ~ 0.02. 

The preparation procedure described in was fol- 
lowed. The liquid crystal was added to the dried particles 
at room temperature, i.e. while in the nematic phase (one 
does not expect good mixing in this situation). We then 
raise the temperature to well above T n j and subject the 
sample to the ultrasonic excitation. In this way the parti- 
cles were gradually dispersed in the isotropic phase of 5CB 
or MBBA; the samples were stored at T ~ 45C (above T n j) 
in a tumbling device to ensure that the mixtures were ho- 
mogeneous before any experiment was started. The par- 
ticle concentration & has been measured by weight. Al- 
though the theoretical arguments require a proper (global 
average) volume fraction <P = Nvr/V, we consider the 
two adequately close because the densities of PMMA and 
of nematic liquid are not very different. Differential scan- 
ning calorimetry (DSC, Perkin-Elmer Pyris 1) was used 
to identify the phase behaviour of the resulting mixtures. 

For this work the major advantage of the confocal opti- 
cal microscope is that it can produce a three-dimensional 
image of a relatively thick and optically inhomogeneous 
sample. Ordinary polarised microscopy, commonly used 
in studies of liquid crystals, does not produce results here 
because the aggregated colloid samples strongly scatter 
light - the samples are opaque. The employment of a pin- 
hole to eliminate the out-of-focus light from planes above 
and below the focal plane in a confocal microscope pro- 
duces a sharp picture of optical contrast in this plane. We 
use a Laser Scanning Confocal Microscope (LSCM 510, 
by K. Zeiss) in the reflection mode which does not require 
fluorescent labelling, illuminated by the monochromatic 
540 nm laser, without polarisers. In this mode, the images 
could be collected at a depth of up to 100 /xm under the 
top sample surface. 



3 Experimental results 
3.1 Calorimetry 

The starting point of practically any theoretical descrip- 
tion of thermotropic nematic colloid phase behaviour, e.g. 
[Q, is the assumption that the nematic transition tem- 
perature in the homogeneously mixed colloid is a linearly 
decreasing function of particle concentration <P. In fact, 
below we shall explicitly calculate the slope of this lin- 
ear dependence: ~ T* (1 — a<P ), eq. (||). Of course, such a 
linear variation means that the particles are regarded as 
fully independent, non-interacting - which cannot be true 
for high concentrations. Also, this estimate for T n ;(^) is 
based on a value reversal of the colloid parameter WR/K 
near the transition. This requires a subtle balance of mate- 
rial constants. If WR/K remains large deep in the nematic 
phase, the particles will quickly ag greg ate under the action 
of strong attraction forces, see Q,[l(j|. If WR/K remains 
small even near T n j in spite of its tendency to grow as the 



4 



V.J. Anderson et at: Cellular solid liquid crystal colloids: 1. Phase separation 




25 30 35 40 25 30 35 40 

Temperature ( C) ^ Temperature ( C) ^ 



Fig. 2. Differential scanning calorimetry of nematic transition 
on cooling a nematic colloid mixture from its homogeneous 
isotropic state. The heat flow peaks are exothermic, cooling 
rate 5°/min. (a) Data for 5CB at concentrations # = 0, 3%, 7% 
and 15% (consecutive curves with decreasing T n i). One should 
note both the shift of the transition temperature and the 
broadening of the peaks with increasing colloid concentration, 
(b) Data for MBBA at concentrations <P = 0, 4.5%, 8% and 
14.5% (also consecutive curves with the peak moving from right 
to left). 
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Fig. 3. The variation of nematic transition temperature with 
particle concentration (average volume fraction $), obtained 
from the DSC analysis. The data for T ni (<£) for 5CB and MBBA 
nematic colloids are fitted with the linear functions: 308.8(1 — 
0.06$) for 5CB - and solid line, and 310.3(1 - 0.16$) for 
MBBA - • and dashed line. 

nematic order parameter Q{T) decreases (cf. section |J), 
the particles will not have a strong elastic energy around 
them and the system will behave more like a molecular 
mixture, a nematic with impurities - a well studied sub- 
ject with a different phase morphology Jl4) . 

In order to determine the coupling energy between the 
nematic field and the individual colloid particle we study 
the isotropic-nematic phase transition on cooling from a 
homogeneously mixed state with a small particle volume 
fraction <£. The differential calorimetry shows an exother- 
mal latent heat peak at this weakly first-order transition. 
Figure || shows the evolution of this peak on increasing the 
particle concentration, in 5CB and MBBA colloids. Two 
main features have to be noted. First, the dependence 
of T n i (determined by the peak onset) on concentration 
shows a monotonic trend, fig. |3[ The linear fits appear 
to be quite good, with the slope consistent with theoret- 
ical estimates of section ||. Secondly, one finds that the 
peaks become broader and of lower total area on increas- 
ing the particle concentration. This is not unexpected, the 
first order transition should become more and more dif- 



fuse with added impurities. In fact, there is a theoretical 
view |l^,[l6| that such a transition should become contin- 
uous in equilibrium. We did not address this question in 
detail, which would have required significantly lower DSC 
cooling rates, at the very least. 



3.2 Cellular structure 




(a) m 




Fig. 4. Reflection-mode confocal images of cellular structures. 
The micrographs show a sample cross-section at the depth ~ 
50 fim, for the 5CB-based colloid at room temperature and 
particle concentrations <P = 3% (a), 7% (b), 10% (c) and 15% 
(d). 



Here we present the confocal images of horizontal cross- 
sections of aggregated nematic colloids. The samples stud- 
ied here were prepared by simply depositing a portion of 
isotropic homogeneous mixture (at 45C) on a slide and 
allowing it to cool to the room temperature, without top 
cover. The contrast mechanism between the nearly pure 
nematic regions and the densely aggregated PMMA par- 
ticles is due to a different average refractive index in the 
reflection mode of the microscope. 

The cellular morphology of aggregated nematic col- 
loids is apparent in fig. ^ for the 5CB system and in fig. || 
for the MBBA system. The relatively large cavities are 
separated by thin interfaces providing an optical contrast. 
The proportion of volume between these two fractions in- 
dicates that the cavities are filled with the nematic liq- 
uid (the majority phase) and the interfaces are made of 
PMMA particles. 

The sizes of cells appear very irregular. One of the 
reasons is that the thickness of confocal imaging plane is 
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Fig. 6. Reflection-mode micrographs of free surface of aggre- 
gated colloids in ambient open-air conditions. Left: 5CB sample 
with <P = 10%, right: MBBA sample with <P = 8%. 



perse, their areas crossing the focal plane would appear 
broadly distributed. Although there is a clear decrease in 
cell size with increasing particle concentration in fig. ^, 
the distribution is too irregular and asymmetric to allow 
accurate quantitative conclusions. 



(c) 



m 



Fig. 5. Reflection-mode confocal images of cellular structures. 
The micrographs show a sample cross-section at the depth ~ 
50 /im, for the MBBA-based colloid at room temperature and 
particle concentrations $ = 2.7% (a), 4.5% (b), 10% (c) and 
14.5% (d). 



much smaller than the cells and one obtains their differ- 
ent cross-sections. However, it is also clear that cells are 
very polydisperse both in size and in shape. One may ar- 
gue that the way of sample preparation has caused the 
polydispersity: the relatively slow cooling under ambient 
conditions passes through the binodal regime of nucleation 
and growth in the phase diagram, see later (fig. ^|), and 
would result in a variety of sizes of growing nematic nu- 
clei. This should be contrasted with a rapid cooling deep 
into the spinodal decomposition regime, where a charac- 
teristic size is selected. We shall see later that samples 
cooled in the rheometer (at a rate of 30°/min, with an 
additional small vibration aiding interface packing) have 
a significantly more regular cellular superstructure. 

Another fact, evident in the confocal images, figs. [| 
and ||, is that the cell walls are perforated. This is espe- 
cially evident in the image (b) for the 7% colloid, where 
most of the walls are continuous in the cross-section, but 
some are clearly interrupted. This observation brings our 
system into the class of "open-cell" structures, cf. fL7| . 
Figure || shows the microstructures of a free surface of ag- 
gregated cellular nematic colloid prepared in the ambient 
conditions. The difference with the internal cellular mor- 
phology is striking. Nevertheless, in spite of a much more 
regular surface, representing an effective cell wall meeting 
the outer interface, the perforations of the interface are 
evident. 

Image analysis, which should provide the statistical 
distribution of cell sizes, is difficult for the confocal scans 
of different cell cross-sections. Even if the cells were monodis- 
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Fig. 7. Analysis of the average cell size A as function of ini- 
tial particle concentration, at room temperature, for PMMA 
particles of R = 150 nm dispersed in (a) 5CB and (b) MBBA. 
Solid lines show the fits to A = 4.88/<5 or n « 32 (5CB), and 
A = 2.22/0 orn w 15 (MBBA). 



4 Theoretical model 

First of all, let us re-iterate that the proposed regime of 
phase separation as a mechanism of particle aggregation 
in the nematic phase is taking place due to the smallness 
of the liquid crystal colloid parameter WR/K. In the op- 
posite case, the study of which has been attempted several 
times over the years, the topological defects around par- 
ticles result in strong interaction forces. If there arc only 
few large particles in the system, they often form "strings" 
- chains of alternating spheres and topological defects (3, 
0| . When the particle concentration is finite, the colloid 
undergoes fast aggregation into 3-dimensional floes. Big 
fiocs with stronger relative anchoring (planar on silica sur- 
face, as in fall into the category of WR/K ^> 1. An 
intermediate stage of evolution is the network of disclina- 
tions connecting the floes and having an effective modulus 
G' ~ 0.01 Pa. The topological defect network would even- 
tually reach the global equilibrium by clearing the nematic 
volume altogether. (This final state is the eventual fate of 
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all other metastable modes of aggregation, including our 
cellular solids). 

We now return to the case of weak anchoring, WR/K <C 
1, and consider the behaviour of a system of small parti- 
cles near the nematic transition, where the amplitude of 
the nematic order parameter Q(T) becomes small (the N- 
I transition is weakly first order, so that the latent heat 
and the jump of Q at the transition are small). Both rel- 
evant elastic parameters, the Frank constant K and the 
anchoring energy W vanish as Q — > 0, but with different 
rates: K ~ kQ 2 , while W ~ wQ in the first approxima- 
tion. Therefore, in the immediate vicinity of T n j the ratio 
WR/K oc 1/Q should become large and the topological 
regime of nematic director distortions around even very 
small particles may prevail. This would certainly be the 
case if a critical second-order transition occurred at T n j. As 
it is, since Q never reaches zero, only a certain combina- 
tion of material parameters would allow WR/K to change 
from small to large. [] The elastic energy of these distor- 
tions is, near the N-I transition, ~ 10KR = IOkRQ 2 . This 
makes an addition to the Landau free energy density de- 
scribing the nematic transition, for local particle concen- 
tration it takes the form ~ 10KR4>/vb,-, with vr 
the particle volume, giving 



F n = ^A o [T-T c (4>)]Q 2 -\BQ* + h0Q i (2) 

Here the shift in the critical point is due to the elastic 
energy around the particles: 



T c (0) =T*[l 



lOnRqh 
A T*v R 



1 
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where the second, very approximate, equation uses the 
fact that the nematic correlation length £ is actually given 
by the ratio of parameters k and A Q , see Appendix. The 
experimental data in fig. || show that the slope of T c (<&) 
is indeed noticeable. The experimental estimate obtained 
from the linear fit is not wholly unreasonable, albeit slightly 
larger than one could expect for our particles and a typical 
correlation length £ ~ 10 nm. 



4.1 Phase separation below T n - S 

Immediately after the transition, the growing nematic nu- 
clei start expelling the particles in the drive to reduce its 
thermodynamic free energy (||). The initial average parti- 
cle volume fraction # splits into locally different valuesin 
phase separating regions; we now use the notation for 
the local concentration. The equilibrium nematic order 
parameter of the mixed phase is 




(4) 



1 If WR/K remains small throughout the nematic phase, the 
elastic energy of deformations, ~ W 2 R B / K ||, has only a weak 
dependence on the nematic order parameter Q and results in 
long aggregation times. 



with the reduced temperature r = TIT* and the coeffi- 
cient a ~ 0.1, using the data of fig. ||. Accordingly, the 
upper limit of stability of homogeneously mixed nematic 
phase is at the concentration 



1 



a \A(A T*)C 



(5) 



The equilibrium free energy density F n (Q*) then becomes 
a function of local particle concentration 0, see eq. (^) 
below. The principal feature of this optimised free energy 
density of a nematic mean field is the rapid increase of 
F n ((j>) with increasing particle concentration and the non- 
convex form of this dependence. This feature is the driving 
force for phase separation. 

This energy penalty on having the particles uniformly 
dispersed in the nematic matrix has to be added to the 
free energy of mixing. The Carnaham-Starling excess free 
energy density of a hard-sphere suspension Jl8j is a very 
good approximation of the equation of state: 



kT 

F p = — [ <j> hxi 
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F p = — 1.790 + 30 In 

Vr 



1 



, , . r < 0*, liquid (6) 
> 0* solid 



where vr is the particle volume and the two expressions 
are matched by an adjustable parameter u « 1.7929 at a 
concentration 0* w 0.52, with C w 0.64 the random close- 
packing fraction. The full free energy of nematic colloid, 
F n + F p , then becomes, in dimensionless form, 



VR 

kT' 
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1 



at small concentrations (0 < 0*), where the dimensionless 
parameters ai 3> 1 and 62 "C 1 are estimated in Appendix, 
eq. (fll|). The factor (1 — 0) in front of the nematic mean 
field energy F n (Q*) measures the local proportion of the 
solvent in the colloid mixture. 

This is the simplest possible model of nematic col- 
loid phase ordering, essentially the same as used in [jjj, 
with minor deviations in the interpretation of nematic 
order-composition coupling and material constants. This 
interpretation, however, determines the results and pre- 
dictions. All parameters entering the free energy F{<j>) will 
be determined from experiment, thus providing the phase 
diagram in the temperature-composition plane with no 
adjustable free parameters. 

At low initial average particle volume fraction, the 
structure of eq. (]?]) is reminiscent of the Flory-Huggins 
model for demixing. The nematic mean-field energy plays 
the role of the effective x-parameter potential term in that 
model, providing non-convex variation in F((f>). Figure |^ 
shows a series of plots of eq. (0) for the set of Landau 
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Fig. 8. Plots of free energy density (^) for decreasing temper- 
ature: T > T* « 309 (Q = 0, thin solid line, eq. |), T = 305 
and 300°K. The nematic energy F n in (Q) has been scaled by 
a factor ~ 10~ 3 to let it fit on the same graph. At concen- 
tration <j>2(T) the nematic nematic transition takes place and 
the isotropic state with Q = takes over. The inset shows the 
fine structure of the plots very near <j> — with a well defined 
minimum at <f>o. 

parameters corresponding to 5CB nematic liquid crystal, 
eqs. ([l0]) and ( |TT| ) in Appendix. There is always a well- 
defined minimum at very small concentrations, 



t>Q ~ exp 



-a 2 



(1 - r)(l - r + 2a) 



due to the entropy contribution (cj> In cj>) in eq. (|6j). The 
usual common-tangent construction for the F(<p) connects 
two branches of the free energy: one non-convex, at very 
low <fr, in an almost pure liquid crystal, and the other fully 
convex, in the high-</> phase where the nematic order is 
unstable (Q = 0) and the particles are densely compacted. 

The schematic phase diagram of the nematic colloid is 
shown in fig. ^ in coordinates of reduced temperature r 
and concentration <j>. At high temperature the colloid be- 
haves as a standard athermal hard-sphere system, which 
crystallises above <fi ~ 0.5 and reaches the maximum ran- 
dom close-packing level at 4> c ~ 0.64. At an extremely low 
particle concentration, corresponding to the expanded in- 
sert in fig. ||, there is a region of homogeneous nematic 
mixture, labelled NL in the phase diagram. This region is 
artificially expanded in fig. the real numbers make it in- 
visible in the plot. A very good analytical approximation 
for this phase boundary is 



Mr) 



1 



a 2 a 



10" 6 t.. 



(8) 



The spinodal on the high-</> branch is determined by 
the point of contact between the nematic (concave, locally 
unstable) and the isotropic (fully convex) parts of the free 
energy, which is marked as 4>2{t) in fig. g. Its analytical 
estimate, 

1 /2, 

2 = - ~& 2 + 1 - T 

a \ 9 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

concentration <|) 

Fig. 9. The phase diagram of a nematic colloid, for the choice 
of parameters discussed in Appendix (thin solid lines) and for 
an artificially small nematic mean field, scaled by ~ 10~ 3 to 
make the phase boundaries apparent (solid binodal lines). In 
the isotropic phase, the hard-sphere colloid is in the liquid state 
(IL) below the coexistence boundary <j> ~ 0.5 and is in the solid 
state (IS) at high cj> — > <j> c . See text for explanations of other 
phase boundaries. 



is practically indistinguishable from the stability limit 
in eq. (|J). This expresses the known fact about the "weak- 
ness" of first order nematic transition: the width of ne- 
matic coexistence is very small in ordinary thermotropic 
liquid crystals. Below the spinodal line 4>2{t) the homo- 
geneous particle mixture is unstable. Above this line, the 
shaded region in fig. || shows the zone of coexistence be- 
tween the \ow-(f> (nematic) and the high-</> (isotropic) phases. 
The upper boundary of this region, the high-concentration 
binodal line is calculated numerically from the common- 
tangent condition for F{4>) in fig. |§|, for the set of param- 
eters given by eq. ( |TT| ) (thin solid line) and the artificially 
scaled down nematic energy (bold solid line) to expose the 
critical point, where this binodal merges with the low-(/) 
phase boundary 4>\{t). 

Note that the high-concentration phase is in the state 
of colloidal solid. This is a consequence of very large rela- 
tive strength of nematic mean field, expressed by the large 
dimensionless parameter in eq. (pT|). One may say that the 
effective pressure from the nematic liquid to expel the par- 
ticles and compact them in high-0 regions is rather high. 

A word of caution has to be sounded at this point, re- 
garding the continuing usage of the linear law T c (<f)) at rel- 
atively high particle concentrations. This approximation 
is not likely to be valid in the regime of high interparticle 
interactions and, therefore, the spinodal 4>2(i~) would de- 
flect down from the straight line in fig. |j. We, however, are 
mostly interested in small particle concentrations where 
the approximation of independent particles would hold. 



4.2 Cellular structure 

We thus envisage the following mechanism of phase evolu- 
tion of the nematic colloid on its cooling from the isotropic 
homogeneously mixed phase: 
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Fig. 10. Model of particle phase separation showing the homo- 
geneous mixing in isotropic phase, the growing nematic nuclei 
expelling the particles and the final state, where the particles 
are densely packed on thin interfaces under the effective pres- 
sure from the nematic mean field in the encapsulated domains 
of size A. 

— After quenching the homogeneous colloid suspension 
below the nematic transition point, the system phase 
separates into growing regions of pure nematic liquid, 
from which the particles are expelled into the bound- 
aries. 

— On these interfaces of growing nematic nuclei, the con- 
centration of particles is so high that the colloid so- 
lidifies and the remaining mesogenic liquid is in the 
isotropic state. 

Clearly, the fine balance of kinetic effects is essential: the 
phase separation should have a higher rate to develop 
within growing nematic nuclei. The driving force for this 
process is the thermodynamic free energy gain in hav- 
ing the high nematic order in a clean system at a given 
T < T n j, in comparison with the uniformly dispersed ne- 
matic colloid, which is expressed by non-convexity of plots 
in fig. |. 

Eventually, when the clean nematic regions grow to 
come into the near-contact with each other, the particles 
are only allowed to pack on narrow interfaces. The sam- 
ple transforms into a cellular superstructure with rather 
thin densely-packed walls encapsulating nematic volumes 
of characteristic size A. This cell size is determined by the 
initial average colloid concentration <P and the thickness of 
interfaces. Let us define this thickness as d = nR, several 
times the particle radius (the smallest possible is n ~ 2 in 
fig. 10). Then the total area of such interface in the sam- 
ple volume V is A = <PV/(nR). This leads to the order of 
magnitude estimate for the mesh size 



(9) 



This gives A ~ 6 /jm for this perhaps unrealistic case of 
n = 2, and R — 150 nm and <P = 0.05. The analysis of 
structures in the next section obtains n ~ 20, see fig. ^, 
i.e. interfaces of ~ 10 particles thick and A ~ 60 fim for 
<P = 0.05. 

The cellular solid with thin high-tension interfaces is 
not the globally equilibrium state. It is obvious that the 
free energy can be further lowered by aggregating all par- 
ticles into a 3-dimensional densely packed group, leav- 
ing the whole volume for the pure nematic liquid phase. 
However, the metastable cellular structure may become 
frozen, trapped by high energy barriers for disrupting the 
formed interfaces. Several mechanisms contribute to creat- 
ing such barriers. Possibly the largest contribution arises 



Fig. 11. An illustration of topological barrier for the breaking 
of interface between two cells, sufficiently big to sustain the 
topological charge, A ^> K/W , see text. 



from the nematic topological argument similar to that in 
liquid crystal emulsions [|19f . Consider the nematic align- 
ment in one cell, regarding the interface as a continuous 
wall imposing a homeotropic director anchoring energy 
W . When the mesh size A increases so that the familiar 
dimcnsionless parameter (WX/K) ^> 1, the nematic liq- 
uid crystal within this cell has to possess the topological 
charge N = 1. When one attempts to break through the 
cell interface, two topologically charged volumes come into 
contact and their coalescence must comply with the law of 
this charge conservation. Therefore, in the moment of for- 
mation of a channel between the cells, a new topological 
defect of the charge N = — 1 must be created [[19| , Later 
this defect may travel towards one of existing monopoles 
and annihilate it to minimise the elastic energy of the fi- 
nal joint volume. However, this cannot happen instantly 
and the formation of each new topological defect costs 
an elastic energy ~ KX. This energy barrier can be very 
high: ~ 10~ 16 J, compared with k^T ~4x 10 _21 J, and 
the system may remain trapped in the metastable state 
of random cellular structure. One then expects to find the 
effects of ageing and time translation invariance breaking, 
characteristic of weakly non-ergodic glassy dynamics and 
rheology MM. 



5 Conclusions 

We have reported the results of a structural study of ne- 
matic colloids based on a classical thermotropic liquid 
crystals matrix with small monodisperse polymer parti- 
cles suspended in it [jyj. The hydrophobic sterical sta- 
bilisation ensures the radial anchoring of the director on 
particle surface. At relatively small concentrations of par- 
ticles, we observe good mixing above the clearing point 
of isotropic-nematic phase transition T n i and a rapid ag- 
gregation of the homogeneous mixture into a rigid gel-like 
solid, completely opaque optically at and below T n {. The 
properties of the phase ordering and the morphology of 
the phase-separated aggregates were the primary focus of 
this study. 

On cooling from the homogeneous isotropic mixture, 
we observe a decrease in the transition temperature T n j as 
a function of average colloid concentration <£, which fol- 
lows a reasonably linear law in the region of small concen- 
trations studied here. This is an expected and frequently 
observed phenomenon, accompanied by a noticeable dif- 
fusion of the weak first-order phase transition, which is 
otherwise sharp in a pure nematic liquid crystal. 
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Below T n ; the colloids undergo phase separation with 
the resulting structure best described as an open cellular 
solid, with the particles densely packed in thin walls and 
the cavities filled with a pure liquid crystal. The character- 
istic cell size is of the order of 10 — 100 /im and appears to 
follow the inverse-proportionality law A ~ l/#, also con- 
firmed by a simple theoretical estimate. Since the nematic 
director is anchored on cell walls, the resulting randomly 
quenched birefringent texture strongly scatters light, giv- 
ing the material its opaque appearance. The remarkable 
mechanical rigidity of the resulting cellular solids, stem- 
ming from the high effective surface tension of interfaces, 
is the subject of the companion paper [^2[. 



Colloid suspensions in nematic liquid crystal matri- 
ces, thermotropic and lyotropic, have been studied before. 
Particle aggregation is an unavoidable result in all cases, 
since the system below T n ; tends to minimise the elastic 
energy of director distortions around individual particles. 
However, the cellular solid morphology, first observed in 
PJ , is new. We believe that two main factors contribute 
to the reason why we obtain such an effect. The parti- 
cles have to be small and have sufficiently low anchor- 
ing energy to ensure that the colloid parameter WR/K is 
small just below T n j. On the other hand, our theoretical 
model suggests that this parameter has to become sub- 
stantial again, to provide a required thermodynamic force 
for phase separation. This rather restricts a range of par- 
ticle sizes and surface treatments that allow for such a 
compromise to occur. Secondly, we find that the rate of 
cooling through the first-order nematic transition has an 
important effect. The cellular structure is much more reg- 
ular and robust when the homogeneous isotropic colloid 
is quenched rapidly deep below T n {. The theoretical phase 
diagram suggests that in order to achieve a good selection 
of size one needs to avoid a binodal region just below T n j. 
In many cases of slow-rate cooling and a colloid material 
parameters making the binodal gap wider (such as in our 
example of MBBA system), the mixture would separate in 
a different fashion - perhaps totally avoiding the cellular 
solid regime. 



Although the details of aggregation mechanism and 
forces that hold particles together at interfaces are some- 
what unclear, the suggested theoretical "toy model" offers 
a possible explanation for the formation of cellular struc- 
ture. Still, many questions remain open. Such an unusual 
phase behaviour and remarkable rheological properties of 
the liquid crystal colloid suspensions require further de- 
tailed study, both theoretical and experimental. 



We appreciate valuable discussions with M.E. Cates, S.M. 
Clarke, P.D. Olmsted and M. Warner. The help of I. Hop- 
kinson with the confocal microscope is gratefully acknowl- 
edged. This research has been supported by the UK EP- 
SRC. 



Appendix: 

Parameters of Landau free energy 

It is useful to find the values of the three phenomeno- 
logical parameters describing the transition, coefficients 
A , B and C in eq. @, although different materials will 
have these parameters slightly different. It is nevertheless 
instructive to examine the characteristic orders of magni- 
tude. In order to determine three parameters one needs 
three independent measurements. We are fortunate that 
there are, in fact, four available: 

1) The jump of order parameter at the weak first-order 
transition is AQ n i = 2B/3C. There might be some error 
in its determination, which should depend (among other 
factors) on the rate of cooling through the transition. How- 
ever, because Q only varies between and 1 and because a 
number of molecular theories predict this jump explicitly, 
one can take qualitatively Q n ; ~ 0.4, giving B » 0.6C. 

2) The second measurement can be the interval be- 
tween the transition temperature and the critical point 
T* (the latter may be determined by extrapolating the 
inverse susceptibility, ~ \T — T*\). In usual ther- 
motropic nematic liquid crystals this interval is rather 
small, T ni — T* = 2B 2 /9A C ~ 1°. This gives B w 
7.5A Q ■ 1°K. (Note that the dimensionality of B and C is 
energy density, while A a has the dimensions of J/m °K). 

3) The latent heat of the first-order phase transition 
is AH = T ni (2A Q B 2 /9C 2 ), per unit volume. It is typi- 
cally obtained from calorimetry by integrating the char- 
acteristic peak. Again, a large uncertainty may accom- 
pany such a measurement because at any non-zero cooling 
rate the conditions are not exactly equilibrium. However, 
keeping up the qualitative approach, the value of transi- 
tion enthalpy per unit mass is of the order ~ U/g for 
5CB at a reasonably slow cooling rate. Taking the density 
p ~ 1 g/cm 3 and the T n ; ~ 310°K, we obtain 



i G ~6x 10 3 J/m 3o K 
Then B~5x 10 4 and C ~ 1.2 x 10 5 J/m 3 



(10) 



4) It is important that the three measurements above 
give estimates that agree with the fourth way of accessing 
A . The nematic correlation length £ may be determined 
by the ratio of the bare Frank elastic constant k = K/Q 2 
to the thermodynamic energy density: £ 2 = k/A d AT. 
There are many ways of confirming that the characteristic 
magnitude of £ is ~ 10 nm. For usual values k ~ 10 _11 J/m 
and AT ~ 10° we obtain A Q ~ 10 4 J/m 3 °K, m a very rea- 
sonable agreement with the previous estimate. 

Finally, we may estimate the dimensionless factors en- 
tering the final expression for the free energy (^) . At T* » 
309°K (kT* ~ 4.7 x 10~ 21 J) and particle volume v R ~ 
1.4 x 10~ 20 m 3 we obtain 



B 2 



«2 



- 2 x 107 and 62 a t r 



10" 



(11) 
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